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Day to day to electrical power demand increases very rapidly with linear and 
non-linear load demands. Especially the nonlinear loads are creating the 
harmonics in the current and voltage signals. The current and voltage signal 
values are measured with the phasor measurement unit (PMU) for the proper 
magnitude and phase angle calculation even in the presence of harmonic 


components in the signals. The performance of the PMU is depending upon 


the phasor calculation technique. Different technique/methods are available 
Keywords: for the phasor calculation, from the method to method there is difference in 
the accuracy, phasor computation time and complexity. In this paper various 


Direct phasor calculation techniques for phasor calculation are presented. For better performance of 


Full cycle DFT PMU, more accurate and less computation time for phasor calculation 
Global positioning system technique is required. But in real time, accuracy and speed both may not 
Half cycle DFT satisfied with single technique. Need to find a satisfactory technique, which 


Least square estimation 
Phasor measurement unit 


satisfies the speed of phasor computation and accuracy. In this paper it is 
proposed that direct phasor estimation technique, which gives the better 
results in terms of accuracy and time and this method, satisfies the 
requirements for the dynamic monitoring of power system according to the 
IEEE std. C37.118.1-2011. 
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1. INTRODUCTION 

For the real time power system monitoring, phasor calculation is very important [1], [2]. Because 
when the electrical power generation and load demand mismatch then there is a change in the magnitude of 
the voltage phase angle and frequency. For satisfactory operation of power system all the above variables 
should monitor properly [3], [4]. For this, phasor measurement unit (PMU) is a perfect device to measure all 
the above/mentioned parameters [5], [6]. The performance of the PMU is evaluated with the phasor 
calculation technique [7], [8]. In real time different methods are there for the phasor calculation. They are 
discrete fourier transform (DFT), least square estimation, cosine transform and zero crossing detection [9], 
[10]. The DFT further classified into full cycle DFT and half cycle DFT based on the signal used for the 
phasor computation [11], [12]. The DFT technique is also classified into two types based on the phasor 
computation process; they are recursive DFT (RDFT) and non recursive DFT (NRDFT) [13], [14]. RDFT 
takes less computation time for the phasor estimation [15], [16] but if any error is occurred in the calculation, 
then it is carried forward to the next calculations [17], [18]. The NRDFT gives the accurate values than 
RDFT [19], [20] but it takes more time for the computation [21], [22]. According to IEEE C37.118.1-2011 
for the dynamic state monitoring of the power system more time for data computation is not suitable [23], 
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[24]. To overcome above problems direct phasor estimation technique is proposed in this paper. This method 
produces better results and total vector error (TVE) also less as compared with other methods. 


2. POWER SYSTEM PHASOR ESTIMATION TECHNIQUES 
For the power system monitoring, power flow calculations the correlation among the signals is very 
important. Two real time signals “X(t)” and “Y(t)” correlation is as follows 


Cxy(m) = 2 X(n)Y (n-m) 


n=- 


Here X(n) and Y(n-m) are the real signals in discrete form, ‘m’ represents the lag in the signal. 
Differences in the phase angles of a same frequency signals can be represented in the single complex polar 
plane with phasor. For the phasor calculation many methods are there, each method performance can be 
calculated with the %TVE [25], [26]. 


(PE — p£) + (p? — p?) 


%TV = 
i (pe)? + (PH)? 


Here, P,°, Pi° are the calculated/estimated phasor real and imaginary values. P,’, P;* is the real/true 
phasor real and imaginary values. For the real and imaginary values calculation different phasor estimation 
techniques are available. 


3. DISCRETE FOURIER TRANSFORM 

Phasor is complex form representation of the sinusoidal signal with magnitude and phase angle 
representation. In the real time many harmonics are there in the current signals due to non-linear loads. 
Consider a real time signal. 


X(t)=ximSin (@1t+@1)+x3msin (@3t+@3)+x5 msin (@st+@5) (1) 
Assume the peak amplitude of the signal is 325; the sampling rate is 12 samples per cycle. The 


sampling of a signal is shown in the Figure 1. By using DFT technique phasor value of a signal can be 
measured. For this half cycle DFT or full cycle DFT can be used. 


Voltage (V) 


0.195 0.2 0.205 021 0215 022 0.225 0.23 
Time 


Figure 1. Half and full cycle window for DFT 


a. Full cycle DFT 
Signal: v(t) = V, sin(wt, + 0) (2) 


Sampling: Vn = VW sin(wt, + 0) 


2mnk 


Voltage phasor, V = ye (vne T) ;0 <n <M- 1, k=0,1,2,3. 
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Where, M=number of samples in a cycle. 


Vn =n" sample of v(t) 


Where k is harmonic component number, k=0 for DC component, k=1 for fundamental, k=2 for second order 
component. The real and imaginary components of the fundamental (k=1) component is: 


V2 aM 2 

Vreat = M neo. (Un cos =) (3) 
v2 œM- . 2 

Vimag = T nao (Un sin =) (4) 


Computed phasor: V = V,eq — jVimag=lVI< 9 


Were, |V] = |Veeat” + Vinag’ 0 = —tan™* (Zmes) 


Vreal 


From (2) Vn = 325 sin(100t, + 30°), sampling rate of 600 Hz, M=12. Full-cycle DFT computation for 
window 1: 


(0.1s to 0.1183s, M=12 points) (0 <n < M — 1). 
V=V2/12[957.224-j1698.8]=229.8002-60.600. 
For window 2: full-cycle DFT computation for window2 (0.1016s to 0.1112s, M=12 points), (O<n<M-1. 
V =a [1678.374-992.588]=229.799 2-30.600. 


The phasor calculation values of both windows are presented in the Table 1. 


Table 1. Full cycle DFT technique phasor calculation values 


Ist window 2nd window 
Time Vv Cos Sin Vin * Cos Vn* Sin y Cos Sin Vn* Cos Vn * Sin 
a (2mn/M)  (2nn/M) (2mn/M) (27n/M) m (27n/M) (27n/M) (2mn/M) (27mn/M) 
0.1 162.5 1 0 162.5 0 - - - - - 
0.1016 277.99 v3/2 %2 240.746 138.99 277.99 1 0 277.99 0 
0.1033 324.98 ” 3/2 162.49 281.440 324.98 V¥3/2 1/2 281.440 162.49 
0.1050 281.45 0 1 0 281.458 281.45 %2 v3/2 140.729 243.749 
0.1066 168.35 -1/2 vV3/2 -84.179 145.802 168.35 0 1 0 168.358 
0.1083 3.403 —v3/2 %2 -2.947 1.701 3.403 -1/2 v3/2 -1.701 2.947 
0.11 -162.5 -1 0 162.5 0 -162.5 —v3/2 1/2 140.729 -81.25 
0.1116 - —v3/2 -1/2 240.746 138.99 - -1 0 277.99 0 
277.99 277.99 
0.1133 - -1/2 —v3/2 162.49 281.440 - —v3/2 -1/2 281.440 162.49 
324.98 324.98 
0.1150 - 0 -1 0 281.458 - -1/2 —V3/2 140.729 243.479 
281.45 281.45 
0.1166 - 1/2 —vV3/2 -84.179 145.802 - 0 -1 0 168.358 
168.35 168.35 
0.1183 -3.403 vV3/2 -1/2 -2.947 1.701 -3.403 1/2 —V3/2 -1.701 2.947 
0.1112 - - - - 162.5 3/2 -1/2 140.729 -81.25 
Summation of the complete cycle samples 957.22 1698.782 Summation of the complete 1678.374 992.588 


cycle samples 


b. Half cycle DFT 


Voltage phasor, V = 


Defining: 
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M/2 


we N2 


n=0 


Mo = je 
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M 
x71 
V2 2an 
Vreal > (Ye S an) 
M/2 M 
n=0 
ae 
2 nn 
Vimag M/2 (Yn sin | 
n=0 


Computed phasor: 


ý= Vreal — JVimag=lVIK 0 


Where |V| = Siar + Vineg? and @ = — tan“? (2) 
Vreal 


From (2) Vn = 325 sin(100t, + 30°). For window 1: Half-cycle DFT computation for window] 
(0.1s to 0.1083s, 6 points), M=12. 


v-22 [478.61-j849.31].V|-229.798. 


849.31 
478.61 


@=tan} (— 25) |V|<0=229.798<-60.99. 


For window 2: half-cycle DFT computation for window] (0.1016s to 0.11s, 6points), M=12. 


ya [839.187-j496.|V|=229.798. 


496.294 
839.187 


@=tan71(— ).|V|KO=229.99<-30.99. 


The calculated phasor value using halfcycle DFT is presented in the Table 2. 


Table 2. Halfcycle DFT technique phasor calculation values 


Ist window 2nd window 
Time Vv Cos Sin Vin *Cos Vn *Sin Vv Cos Sin Vin * Cos Vn * Sin 
a (2mn/M) __(2an/M) (2an/M) (27n/M) "k (27n/M) (27n/M) (27n/M) (27n/M) 
0.1 162.5 1 0 162.5 0 
0.1016 277.99 v3/2 1/2 240.746 138.99 277.99 1 0 277.99 0 
0.1033 324.98 %2 v3/2 162.49 281.440 324.98 vV3/2 1/2 281.440 162.49 
0.1050 281.458 0 1 0 281.458 281.458 %2 V¥3/2 140.729 243.749 
0.1066 168.358 -1/2 3/2 -84.179 145.802 168.358 0 1 0 168.358 
0.1083 3.403 —V3/2 1/2 -2.947 1.701 3.403 -1/2 V¥3/2 -1.701 2.947 
0.11 - - - - - -162.5 -¥3/2 %2 140.729 -81.25 
Summation of the complete cycle samples 478.61 849.391 Summation of the complete cycle 839.187 496.294 
samples 


Even in the signal harmonic content is there, but with this phasor measurement technique only 
fundamental component magnitude and phase angle is calculated. But this DFT technique will give error 
value when the signal having a DC component. Generally, the decaying component is present in the signal 
due to faults; under this condition DFT based phasor estimation may not give accurate values. 


4. DIRECT PHASOR CALCULATION TECHNIQUE 

A new approach is proposed to calculate the instantaneous phase angle and magnitude of the voltage 
signal. This method will give instantaneous phase angle and voltage values more accurately with time 
stamping according to IEFE standard C37.118.1-2011 for dynamic monitoring of the power system. The 
calculated values are time stamped with Universal Coordinated Time (UTC), with this it is possible to 
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synchronize the data from the different locations of the power system, and dynamic power flow in the lines is 
possible. The flow chart for the direct phasor calculation is shown in the Figure 2. According to the flow 
chart the phasor calculation are implementaed in the LabVIEW software. The input sampled signal, calculate 
phasor values magnitude and phase angleare shown in the Figure 3. The timereference for the phasor values 
calculation is taken from the global positioning system (GPS) module NEO-6M. The real time voltage signal 
GPS signals are interfaced to LabVIEW software through NI-MyDAQ card. The calculated phasor values are 
tabulated in the Table 3. 


Read the sampled input signal 
V0=Previous value of the signal 
V1=Present value of the signal 


V=Vp/Sart(2) 


#=0+180 


Write to measurement file 


Figure 2. Flowchart of direct phasor calculation technique 


Instantaneous, RMS and Phase angle 


= ais aT 
12:14:51.014 P. 


400 
300 
200 
100-43 
o 
100 
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400 


Figure 3. Magnitude and phase angle measured from direct phasor estimation 


The laboratory-based experimental setup is shown in the Figure 4. Experimental results are obtained 
by using three phase system with 6ampers of load current. The sending end voltage and receiving end 
voltages, instantaneous voltage phase angles are monitored in PC using LabVIEW software. The sensed 
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signal phasor values are calculated with different phasor estimation techniques. The calculated phasor values 
magnetudes are validated with the fluke voltmeter readings. 


Table 3. Phasor values of different phasor calculation techniques 


Direct phasor Recursive phasor Non-recursive 
Time Instantaneous Real Imaginary calculation calculation phasor calculation 
values value values ; Phase : Phase Magnitud Phase 
Magnitude Magnitude 
angle angle e angle 
02-16-2022 0 230 0 230.01 0 229.8097 -75 229.8 0 
12:14:53.052 
02-16-2022 100.5137 151.8727 -172.727 230.01 18 229.8097 -75 229.8 18 
12:14:53.771 
02-16-2022 191.1884 -29.4316 -228.109 230.01 36 229.8097 -75 229.8 36 
12:14:53.923 
02-16-2022 263.1482 -190.741 -128.521 230.01 54 229.8097 -75 229.8 54 
12:14:54.111 
02-16-2022 309.3493 -222.467 58.38069 230.01 72 229.8097 -75 229.8 72 
12:14:54.325 
02-16-2022 325.2691 -99.9965 207.1249 230.01 89.985 229.8097 -75 229.8 90 
12:14:54.472 
02-16-2022 309.3493 86.36702 213.1683 230.01 108 229.8097 -75 229.8 108 
12:14:54.637 
02-16-2022 263.1482 217.1164 75.89786 230.01 126 229.8097 -75 229.8 126 
12:14:54.785 
02-16-2022 191.1884 200.3641 -112.935 230.01 144 229.8097 -75 229.8 144 
12:14:54.931 
02-16-2022 100.5137 47.4913 -225.044 230.01 162 229.8097 -75 229.8 162 
12:14:55.068 
02-16-2022 0 -137.646 -184.265 230.01 180 229.8097 -75 229.8 180 
12:14:55.206 
02-16-2022 -100.514 -229.271 -18.3029 230.01 198 229.8097 -75 229.8 -162 
12:14:55.358 
02-16-2022 -191.188 -165.137 160.0934 230.01 216 229.8097 -75 229.8 -144 
12:14:55.481 
02-16-2022 -263.148 11.1858 229.7278 230.01 234 229.8097 -75 229.8 -126 
12:14:55.616 
02-16-2022 -309.349 179.9096 143.2918 230.01 252 229.8097 -75 229.8 -108 
12:14:55.770 
02-16-2022 -325.269 225.7825 -43.8435 230.01 269.98 229.8097 -75 229.8 -90 
12:14:55.892 
02-16-2022 -309.349 119.0932 -196.766 230.01 288 229.8097 -75 229.8 -72 
12:14:56.006 
02-16-2022 -263.148 -69.1297 -219.365 230.01 306 229.8097 -75 229.8 -54 
12:14:56.134 
02-16-2022 -191.188 -210.388 -92.9354 230.01 324 229.8097 -75 229.8 -36 
12:14:56.282 
02-16-2022 -100.514 -208.716 96.63179 230.01 342 229.8097 -75 229.8 -18 
12:14:56.421 
%TVE 0.93 4.06 4.16 


M-G Set 
For Controlled ee a 
Voltage Sours 3-Phase Resistive Load 


Interfacing Circuit Monitoring and Simulation Circuit 


Figure 4. Experimental setup for the different phasor calculation methods 
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5. CONCLUSION 

In this paper it is presented that simple dynamic state power systems direct phasor estimation 
technique. This proposed method is effective as compared with the other methods, according to IEEE Std. 
C37.118.1-2011. This method is allowing the calculation of frequency, magnitude, and phase angle of the 
real time power system signals. The working of the proposed algorithm is tested under various conditions 
like simulate signal in LabVIEW, simulate signal with noise and real-time voltage signal. With the time 
stamping of measured phasor values dynamic power flow studies also very easy and comparatively this 
method gives better results than Recursive DFT and Non-Recursive DFT phasor calculation methods. 
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